Gain-Driven Discrete Breathers in VT— Symmetric Nonlinear Metamaterials 
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We introduce a one dimensional parity-time ("PT)-symmetric nonlinear magnetic metamaterial 
consisted of split ring dimers having both gain and loss. When nonlinearity is absent we find a 
transition between an exact to a broken PT-phase; in the former the system features a two band 
gapped spectrum with shape determined by the gain and loss coefficients as well as the inter-unit 
coupling. In the presence of nonlinearity we show numerically that as a result of the gain/dissipation 
matching a novel type of long-lived stable discrete breathers can form below the lower branch of 
the band with no attenuation. In these localized modes the energy is almost equally partitioned 
between two adjacent split rings on the one with gain and the other one with loss. 

PACS numbers: 63.20.Pw, 11.30.Er, 41.20.-q, 78.67.Pt 



Considerable research effort has recently focused in the 
investigation and developement of artificial materials ap- 
propriately engineered to exhibit properties not found in 
nature. In the electromagnetic domain, these advances 
resulted in the construction of metamaterials, a class of 
novel artificial structures that provide full access to all 
four quadrants of the real permittivity - permeability 
plane [l[. Recently, there has been increasing interest 
in synthetic materials that exhibit a combined parity - 
time (VT) symmetry Although quantum systems de- 
scribed by VT— symmetric Hamiltonians have been stud- 
ied for many years , the notion of VT symmetry has 
now been extended to dynamical lattices, particularly in 
optics, where VT related concepts can be realized and 
tested. Soon after the developement of the theory of 
VT— symmetric optical lattices [3|,[H , the VT— symmetry 
breaking was experimentally observed 0^3- Naturally, 
such considerations have also been extended to nonlinear 
lattices where the existence of stable discrete solitons Q 
and Talbot effects Q was theoretically demonstrated. 

Among many recent developements in VT— symmetric 
materials, the application of these ideas in electronic 
circuits (lfj, not only provides a platform for testing 
the new ideas within the framework of easily accessi- 
ble experimental configurations, but also provides a link 
to the other important class of artifial structures, viz. 
metamaterials. Magnetic metamaterials, in particular, 
exhibit high-frequency magnetism and negative perme- 
ability in certain frequency intervals 11 1. Conventional 
metamaterials comprising metallic elements suffer from 
high losses that hamper any possibility for their use in 
practical applications. However, building metamaterials 
with VT symmetry, relying on gain and loss, may pro- 
vide a way out and moreover lead to new extraordinary 
properties. It is shown that VT— symmetric metamate- 
rials undergo spontaneous symmetry breaking from the 
exact VT phase (real eigenfrequencies) to the broken VT 
phase (at least a pair of complex eigenfrequencies) , with 



variation of the gain/loss coefficient. In the presence of 
nonlinearity, the generation of long-lived nonlinear exci- 



12] it is 



tations in the form of discrete breathers (DBs) 
demonstrated numerically in discrete metamaterial mod- 
els. These novel gain-driven DBs, that may result either 
by proper intialization of the VT metamaterial or purely 
dynamically through external driving, result from power 
matching of the input power through the gain mechanism 
and internal loss. 
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FIG. 1: (Color online) Schematic of a one-dimensional 
VT— symmetric metamaterial. Upper panel: all the SRRs 
are equidistant. Lower panel: the separation between SRRs 
is modulated according to a binary pattern (VT dimer chain) . 



Consider a one-dimensional array of periodically ar- 
ranged dimers, each comprising two split-ring resonators 
(SRRs); one with loss and the other with equal amount of 
gain (Fig. 1). The SRRs are coupled magnetically and/or 
electrically through dipole-dipolc forces [l3|, EH , and are 
regarded as RLC electrical circuits, featuring a resistance 
R, an inductance L, and a capacitance C. The SRRs may 
become nonlinear by the insertion of a varactor diode in 
their slits [lH . In order to achieve both nonlinearity and 
gain, a tunnel (Esaki) diode [3| that features a nega- 
tive resistance part in its current-voltage characteristics 
can be employed [ItJ. The incorporation of active con- 
stituents in metamaterials that provide gain through ex- 
ternal energy sourses has been reco gniz ed as a promising 
technique for compensating losses [l8[ . Beyond merely 
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loss compensation, stable metamaterials with net gain 
have been demonstrated, both as left-handed transmis- 
sion lines [TtJ and optical fishnet structures . 

Within the framework of the equivalent circuit model 
[l4 , 13" 22 1, extended for the VT— dimer chain, the dy- 
namics of the charge q n accumulated in the capacitor of 
the n— th SRR is governed by 



Ajvf?2n + <?2n+l + Am<Z2u+2 + X' E q2n + l2n+l + ^El2n+2 



e sin(fir) - aql n+1 - (3q : 



2n+l - 792n+l 



(1) 



^M<hn-l + <hn + X' M q2 n +l + A_Eg2n-l + <?2n + X' E q 2n +l 

= e sin(^T) - aq\ n - fiq\ n + jq 2n 



(2) 



m, a m and Xe, X' e are the magnetic and electric 



where A m , X' { 

interaction coefficients, respectively, a and /3 are dimen- 
sionless nonlinear coefficients, 7 is the gain/loss coeffi- 
cient (7 > 0), eo is the amplitude of the external driving 
voltage, while 51 and r arc the driving frequency and tem- 
poral variable, respectively, normalized to the inductive- 
capacitive (LC) resonance frequency ljq and inverse LC 
resonance frequency u)q , respectively, ujq = 1/\/LCq 
with Co being the linear capacitance. By substituting 
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„+i = B e *((an+i)«-nr) j (3) 



where k is the normalized wavevector, into Eqs. ([TJ and 
@, and requesting nontrivial solutions for the resulting 
stationary problem, we obtain 



Q 2 = (-6 ± a/A) /(2a) 



(4) 



where 



— 1 — (Am — A' 



Mi 



b = 



2 [1 - (A B - A S )(A M - A' M )] + + e' K fi K , c = 1 - 
(Ab — X' E ) 2 — e K e' K , A = b 2 — 4ac, and e K = 2Xe cos(/c), 
e' k = 2A' b cos(k), ^ k = 2A M cos(/c), // K = 2A' m cos(k). 
In the following, we assume that the electric coupling is 
negligible, i.e., that Xe = X' E = 0, for simplicity. Then, 
Eq. dU reduces to 
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2(1 - (Am - A' M ) 2 - /i KM 'J 
The condition for having real f2 K for any n then reads 
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under the assumption AmA' m > 0. It is easy to see that 
for Xm = X' M the earlier condition cannot be satisfied for 
all re's for any positive value of the gain/loss coefficient 
7, implying that a large VT— symmetric SRR array (Fig. 
1, upper panel) will be in the broken phase. To the con- 
trary, for Am 7^ X' M , i.e., for a periodic dimer chain (Fig. 
1, lower panel), the above condition is satisfied for all k's 
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FIG. 2: (Color online) (a) Frequency bands for a 
VT— symmetric metamaterial with balanced gain and loss 
(7 = 0.002) for Am = -0.17 and X' M = -0.10 (black and 
red solid lines); Am = —0.12 and A^ = —0.10 and (green 
and blue dashed lines). The imaginary parts are zero (not 
shown), (b) & (c) Real and imaginary parts, respectively, 
of the frequency bands for a VT— symmetric metamaterial 
with balanced gain and loss (7 = 0.002) for Am = —0.17 and 
A' AJ = —0.1699. (d) The boundary between the exact and the 



broken phase on the 7 
Am = -0.17. 



vm plane, where ru — X' m /Xm, for 



for 7 < 7c — I Am — X' M \, (j 4 ~ 0). In the exact phase 
(7 < 7 C ), the VT— symmetric dimer array has a gapped 
spectrum with two frequency bands (Fig 2a). The width 
of the gap separating the bands decreases with decreas- 
ing I Am — X' M \ for constant 7. For 7 ~ 7 C the gap closes 
(Fig. 2b), some frequencies in the spectrum acquire an 
imaginary part as shown in Fig. 2c, and the VT meta- 
material goes into the broken phase. The two phases are 
separated by the exceptional point 7 = 7c, which varies 
with the other parameters. A typical phase diagram on 
the 7 — tm plane, with tm = X' m /Xm is shown in Fig. 
2d. For fixed Xm,X' m , the bandwidths as a function of 7 
are shown in Fig. 3. 

Eqs. (fT]) and ^ implemented by qo{r) = qN+i( T ) = 
that accounts for the termination of the structure in finite 
systems, are integrated numerically with initial condition 

q m (0) = (-l) m_1 sech(m/2), q m (0) = 0, (7) 
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FIG. 3: (Color online) Frequency bands boundaries as a func- 
tion of the gain/loss parameter 7 for a VT— symmetric meta- 
material with balanced gain and loss in the dimer array con- 
figuration, for Am = —0.14, and X' M = —0.04. 




£0 = 0, and nonlinear coefficients a — —0.4, f3 = 0.08. 
Using a gain/loss function of the form 



m = 1 , . . . , Nt 

m = Ne + l,...,N-Nj, (8) 
m = N-Nt + l N, 



we obtain gain-driven DBs that remain stable for very 
long times (> 10 s time units) for relatively wide parame- 
ter intervals. Note that the actual VT— symmetric dimer 
chain is of length N ~ 2Ni, while both its ends are joined 
to lossy dimer chains of length Ng. We found empirically 
that this is the most effective way to stabilize DBs in 
the VT metamaterial, i.e., by embedding it into a lossy 
metamaterial. The lossy parts help the excess energy to 
go smothly away during the long transient phase of inte- 
gration, and thus prevents the blowing up of the solution 
that otherwise may have occured. Fig. 4 shows typical 
examples of the energy density evolution of a breather 
excitation in the n — r plane. The largest part of the 
total energy is concentrated into two neighboring sites 
belonging to the same gain/loss dimer. The correspond- 
ing instantaneous current profiles i n as a function of n at 
maximum current are shown in Fig. 5. Note that these 
profiles are neither symmetric or antisymmetric at the 
SRR level. 

The breather oscillation frequency fis can be obtained 
from the power spectrum of a time series of the energy in 
one of the breather sites. Such a typical power spectrum 
is shown in Fig. 6a, where the logarithm of the power 
spectrum is plotted as a function of frequency; strong 
harmonics to the fundamental frequency Qb are clearly 
observed. The frequency VIb increases with decreasing 
I Am — A' M | as can be seen in Fig. 6b, while the variation 
of the gain/loss coefficient 7 within 10 -2 — 10~ 3 has ap- 
parently no effect on fis at least in the particular case 
shown in Fig. 6c. The dependence of on |Aj\/ — X' M \ 
is illustrated in Fig. 6d, for those values of \Xm — A' M | for 
which DB excitations are stable (or at least long-lived). 

For a gapped linear spectrum, large amplitude linear 
modes become unstable in the presence of driving and 
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FIG. 4: (Color online) Spatiotemporal evolution of the energy 
density E n for approximately two LC— periods for N = 70, 
N e = 10, \' M = -0.10, Xe = X' E = 0, 7 = 0.002, and (upper 
panel) Am = —0.17; (lower panel) Xm = —0.21. Energy 
localization at two neighboring sites, one with gain and one 
with loss, is clearly observed. Only part of the array is shown 
for clarity. 




FIG. 5: (Color online) Gain-driven, current breather instan- 
taneous profiles i n as a function of n at maximum current 
amplitude, for N = 70, Ne = 10, X' M = -0.10, Xe = X' E = 0, 
7 = 0.002, and (a) Am = -0.17; (b) Am = -0.21. Only part 
of the array is shown for clarity. 



nonlinearity. If the curvature of the dispersion curve in 
the region of such a mode is positive and the lattice po- 
tential is soft, large amplitude modes become unstable 
with respect to formation of DBs in the gap below the 
linear spectrum [23[ . For the parameters of Fig. 2a, the 
bottom of the lower band is located at f^o — 0.887, where 
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FIG. 6: (Color online) (a) The logarithm of the power spec- 
trum of energy at one of the central breather sites y = 
logu)[PS(E n=N / 2 )] is shown as a function of the frequency 
for JV = 70, X' M = -0.10, Am = -0.13, \ E = X' B = 0, 
7 = 0.002. (b) The same as in (a) around the fundamental 
breather frequency Qb, for Am = —0.21 (blue - long dashed); 
Am = -0.20 (green - dotted); Am = -0.17 (black - solid); 
Am = —0.14 (red - dashed); Am = —0.13 (maroon - dotted 
dashed). The other parameters as in (a), (c) The same as in 
(a) and (b) around the fundamental breather frequency f2s, 
for three different values of the gain/loss coefficient 7: 0.002 
(black); 0.005 (green); 0.01 (red), (d) Dependence of Q,b on 
the magnetic coupling coefficient Am, for the spectra shown 
in (b). 



the curvature is positive. Moreover, the SRRs are sub- 
jected to soft on-site potentials for the selected values of 
the nonlinear coefficients a and (3. Then, DBs can be 
generated spontaneously by a frequency chirped driver 
as it is illustrated in Fig. 7, where the evolution of the 
energy density E n for a T^T— symmetric dimer chain is 
plotted on the n — r plane. We use the following pro- 
cedure: At time r = 0, we start integrating Eqs. ([T]) 
and ([2|) with zero initial state and external driving for 
500 T ~ 3500 time units (t.u.), where T a = 2ir/n , 
to allow for significant developement of large amplitude 
modes. At time r ~ 3500 t.u. (point A on Fig. 7), the 
driver is switched-on with low-amplitude and frequency 
slightly above Ct (1.01 flo — 0.894). The frequency is 
then chirped downwards with time to induce instabil- 
ity for the next 10,600 t.u. (~ 1500 T ), until it is 
well below n (0.997 n ~ 0.882). During that phase, 
a large number of excitations are generated that move 
and strongly interact to each other, eventually merging 




FIG. 7: (Color online) Spatiotemporal evolution of the energy 
density E n for a VT— symmetric dimer chain with N = 70, 
Nt = 10, f2 = 0.887, 7 = 0.002, Am = -0.17, \' M = -0.10 
(A B = \' E = 0), and (a) e = 0.085; (b) e = 0.095. The ver- 
tical lines separate different stages in the chirping procedure 
(see text). 



into a small number of high amplitude (multi-)brcathcrs. 
At time r ~ 14, 100 t.u. (point B on Fig. 7), the driver 
is switched off and the DBs that have formed are solely 
driven by the gain. They continue to interact until they 
reach a stationary state and get trapped at particular 
sites. The high density segments between points B and 
C in Fig. 7 present precisely those stationary gain-driven 
breathers and multibreathers generated through chirping 
and subsequent dynamics. At time r ~ 85150 t.u. (point 
C on Fig. 7), the gain is replaced by equal amount of loss, 
and the breathers die out rapidly. 

The construction of metamatcrial elements with gain is 
nowadays feasible. In the microwave frequency range, in 
particular, the use of negative resistance devices like tun- 
nel diodes [17J, [24[ and transistors [25j seems to be very 
convinient for providing gain in transmission lines and 
SRRs. Thus, a stable VT metamaterial whose dynamics 
is approximately described by model equations like ([T]) 
and j2]) can be realized and tested. In the specific case of 
a balanced configuration, where dissipation is counterbal- 
anced by equal amount of gain, it is demonstrated numer- 
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ically that stable or at least long-lived, gain-driven DBs 
may be generated. DB generation can be achieved cither 
by proper initialization of the system, or purely dynam- 
ically by frequency chirping of a weak alternating driv- 
ing applied homogeneously to all elements of the dimer 
chain. Importantly, single-site DBs apparently cannot be 
formed but instead, fundamental breather excitations oc- 
cupy two neighboring sites, one with gain and the other 
with loss, belonging to the same dimer. Stabilization 
of these localized modes can be achieved by embedding 
the VT— symmetric dimer chain into a lossy dimer chain, 
as it was found empirically. It is also worth mentioning 
that in a slightly unbalanced configuration, where abso- 
lute gain and loss differ by a small amount, DBs may still 
be generated through the chirping dynamical procedure 
that are not stable in the long term. When loss exceeds 
gain, the generated DBs decay slowly until they vanish in 
a time-scale proportional to the gain/loss imbalance. In 
the opposite case, the DBs gain energy in the course of 
time becoming wider while in the long term their energy 
diverges. 
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